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Abstract 

Using data collected with the CLEO III detector at the Cornell Electron Storage Ring, we study 
the inclusive production of baryons/antibaryons (p, A) and mesons (4> and /2(1270)) in gluon- 
fragmentation and quark-fragmentation processes. We first corroborate previous per-event total 
particle yields in T(1S)— > ggg compared with nearby continuum (e + e~ — > qq) indicating greater 
(~ x2) per-event yields of baryons in 3-gluon fragmentation. We find similar results when we 
extend that comparison to include the T(2S) and T(3S) resonances. With higher statistics, we 
now also probe the momentum dependence of these per-event particle yields. Next, we compare 
particle production in the photon-tagged process T(1S) — > ggj with that in e + e~ — > qqj events, 
to allow comparison of two-parton with three-parton particle-specific fragmentation. For each 
particle, we determine the 'enhancement' ratio, defined as the ratio of particle yields per gluon 
fragmentation event compared to quark fragmentation event. Thus defined, an enhancement of 
1.0 implies equal per-event production in gluon and quark fragmentation. In the photon-tagged 
analysis (T(1S) — > ggj compared to e + e~ — ► qq"/), we find almost no enhancement for protons 
(~ 1.2 ± 0.1), but a significant enhancement (~ 1.9 =t 0.3) for A's. This small measured proton 
enhancement rate is supported by a study of baryon production in Xb2 gg P + A relative to 
Xbi ~^ QQ9 p + X. Overall, per-event baryon production in radiative two-gluon fragmentation is 
somewhat smaller than that observed in three-gluon decays of the T(1S). Our results for baryon 
production are inconsistent with the predictions of the JETSET (7.3) fragmentation model. 

PACS numbers: 12.38.Aw, 12.38.Qk, 13.60.Hb, 13.87.Fh 
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I. INTRODUCTION 



Understanding hadronization, the process by which elementary partons (gluons and 
quarks) evolve into mesons and baryons, is complicated by its intrinsically non-perturbative 
nature. Due to the fact that gluons carry two color indices whereas quarks carry only one, 
the intrinsic gluon-gluon coupling strength (Ca=3>) is larger than the intrinsic quark-gluon 
coupling strength (Cjr=4/3). Radiation of secondary and tertiary gluons is therefore ex- 
pected to be more likely when hadronization is initiated by a gluon rather than by a quark. 
This results in a greater number of final state hadrons as well as a larger average transverse 
momentum in the former case compared to the latter case. In the limit Q 2 — > oo, the ratio of 
the number of hadrons produced in gluon-initiated jets to the number of hadrons produced 
in quark-initiated jets is expected, in lowest order, to approach a simple color-counting ratio 
9/40. 

Many experiments have searched for, and found, multiplicity and jet shape differences 
between quark and gluon fragmentation. At Z° energies, qqg events are distinguished by their 
three-jet topology. Within such events, quark and gluon jets can be separated by a variety 
of techniques including vertex tagging. Because gluons rarely fragment into heavy quarks, 
they will produce jets that form a vertex at the e + e~ interaction point. Quark jets, to the 
contrary, tend to form a detached vertex when the jet contains a long-lived bottom or charm 
quark. For light-quark events with gluon radiation, however, the assignment of final state 
hadrons to the initial state partons is generally more ambiguous and often relies on Monte 
Carlo simulations to determine the fraction of times that an observed hadron is correctly 
traced to a primary parton. At lower energies, one can exploit the decay characteristics of 
quarkonium states to directly compare gluon and quark fragmentation using data taken both 
on- resonance and off- resonance (on the continuum), respectively. The 10 GeV center of mass 
energy range offers a unique opportunity to probe quark and gluon fragmentation effects, 
without relying on Monte Carlo simulation to associate the final state hadrons with an initial 
state parton. CLEO[2] found that the thrust and charged multiplicity distributions of Xbo 
and Xb2 two- gluon decays are more similar to T(1S)— > ggg than to continuum e + e~ — > qq 
events; the reverse was found to be true for Xbi ~^ Qqg- 

Specific particle production in gluon- and quark-fragmentation has also been studied. 
Within the limits of their precision, previous studies at SLD found inclusive production of 
pions, kaons and protons to be equivalent for gluon-tagged and quark-tagged jets^. OPAL 
has measured inclusive charm production to be (3.20 ± 0.21 ± 0.38)% in gluon jets0, 0], 
more than an order of magnitude smaller than the rate observed in quark jets at the Z°. 
ALEPH[6] and DELPHI0] both measured inclusive bottom production in gluon-tagged jets 
to be 2 — 3 x 10~ 3 , again considerably smaller than that expected from charge counting 
in quark fragmentation. Most directly comparable to our current work, OPAL has also 
compared inclusive K® and A production in gluon- vs. quark-tagged jets in e + e~ — ► qqg 
events, finding inclusive production ratios (g/q) consistent with unity (0.94 ± 0.07 ± 0.07 
and 1.18 ± 0.01 ± 0.17, respectively) 0] . 

The decay T(1S) — > gg'y allows one to directly compare the gg system in a gg^y event with 
the qq system in e + e _ — > qqy events. In these cases, the system recoiling against the photon 
consists (to lowest order) of hadrons that have evolved from either a two-gluon or a quark- 
antiquark system. The properties of the recoil systems can then be compared. 1 Additionally, 



Although there may be gluon radiation from the initial partons, we do not distinguish such radiation 
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the radiative transitions from the radially excited T states to the orbitally excited Xb triplet 
offer an opportunity to further probe fragmentation differences between decays of the J=0 
and J=2 Xb states, which decay predominantly to two gluons, vs. decays of the J=l state. 
Since the J=l state is prohibited from decaying into two on-shell gluons, the decay into one 
hard and one soft, nearly on-shell virtual gluon (gg*, followed by g* — > qq) is kinematically 
most favored. Statistical correlations between transition photons with inclusive production 
of particular final-state particles (X) allows a measurement of the relative yields of gg — > 
X '■ qq(g) — > X to these species. 

In an over-simplified 'independent fragmentation' model, hadronization occurs indepen- 
dently for each parton. In such a picture, if fragmentation of each parton (gluon or quark) of a 
given energy is identical, then the ratio of particle production for ggj : qqj : (xb —>■ gg) '■ ggg 
hadronization should vary as: 2:2:2:3. In the opposite extreme, fragmentation occurs in the 
stretching 'strings' between the two partons, in which case the above ratio should be 1:1:1:3. 

In this analysis, we focus on the relative production rates of baryons (jp and A — > prt) 
and heavy mesons (<f) — > K + K~~ and / 2 (1270) — > 7r + 7r~) in gluon vs. quark fragmentation 
(charge conjugation is implied). A previous study noted enhancements in the production of 
0, A and p in three-gluon decays of the T(1S)[^], at a statistical significance of no more than 
than 2-3 a. That initial study also found approximately one unit larger charged multiplicity 
for three-gluon fragmentation of the T(1S) compared to qq fragmentation at a comparable 
center-of-mass energy. With the limited statistics available at that time, the additional 
unit of multiplicity could entirely be accounted for by enhanced three-gluonic production of 
baryons. We now have sufficient statistics to re-measure the three-gluon particle production 
rates, and also to compare, for the first time, inclusive production in two-gluon fragmentation 
vs. inclusive production in three-gluon fragmentation. 

Since then, other experimental data on quark/gluon fragmentation differences in the 
yfs ~10 GeV energy regime have become available, including: 

1. The observation that fragmentation of the J=l state of the Xb triplet (xbi — * qqz^ 
charm) results in charm production comparable to the underlying continuum; no such 
charm production is observed in the two-gluon decays of the J=0 or J=2 states 13 ■ 



2. An enhancement in production of hidden charm in gluo nic decays of the T resonances: 
(T(lS)-> ggg -> W + X)/(e + e- -> J/ip + X) >5[ll| at 90% c.l. 

3. Production of deuterons from resonant 3-gluon decays of both the T(1S) and T(2S) 
at the level of 10 -3 ; no significant production of deuterons is observed from the 



continuum [121]. Enhancements per event are >10. 



Production of rj' in gluonic decays of the T resonance of similar magnitude to that 
observed in T decays via qq: {T(1S)— » ggg — > rf + X)/(T — > qq — > rf + X) ~2/3), 
integrated over momentum [13|]. 



explicitly in this analysis. Thus, the states that we are comparing are, strictly speaking, ggj and qq"f 
to lowest-order only; additional gluon radiation, to which we are not experimentally sensitive, may be 
present in many of the events in our sample. Without the ability to adequately identify additional gluons, 
such higher-order radiative effects are therefore implicitly absorbed into the experimental measurement. 
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II. DETECTOR AND DATA SAMPLE 



The CLEO III detector [H, M, M is a general purpose solenoidal magnet spectrometer 
and calorimeter. The main components of the detector used in this analysis are the drift 
chamber and the silicon detector used for track finding, the crystal calorimeter for energy 
measurements, and the Ring Imaging Cherenkov detector (RICH) and specific ionization 
loss in the drift chamber for particle identification. This system is very efficient (e >98%) 
for detecting tracks that have transverse momenta (pr) relative to the beam axis greater 
than 200 MeV/c, and that are contained within the good fiducial volume of the drift cham- 
ber (| cos#| <0.93, with 9 defined as the polar angle relative to the beam axis). Below this 
threshold, the charged particle detection efficiency in the fiducial volume decreases to ap- 
proximately 90% at pr ~100 MeV/c. For px <100 MeV/c, the efficiency decreases roughly 
linearly to zero at a threshold of px ~30 MeV/c. Just within the solenoidal magnet coil is 
the electromagnetic calorimeter, consisting of 7800 thallium doped Csl crystals. The central 
region of the calorimeter covers about three-quarters of the solid angle and has an energy 
resolution of 

with E the shower energy in GeV. This parameterization translates to an energy resolution 
of about 2% at 2 GeV and 1.2% at 5 GeV. Two end-cap regions of the crystal calorimeter 
extend solid angle coverage to about 95% of 47r, although energy resolution is not as good 
as that of the central region. The tracking system, RICH particle identification system and 
calorimeter are all contained within the 1.5 Tesla superconducting coil. Flux return and 
tracking chambers used for muon detection are located immediately outside the coil and in 
the two end-cap regions. 

We use the CLEO-III data collected at the narrow T resonances as a source of ggg 
and gg'y events, and data taken just below the narrow resonances, as well as the below-4S 
continuum (y / i=10.55 GeV) as a source of qq and qq'j events. Since T(4S)— > BB ~100%, 
data collected on the broad T(4S) resonance is analyzed as a 'control' sample, for which 
we expect no deviation from the below-4S continuum when we require a photon having 
z 7 = E 7 /E heam >0.5. 

The 7 in our continuum qq'j sample results primarily from initial state radiation (ISR) [lj 
We compare events for which the fractional photon energies are the same, which ensures 
that the recoil systems (either two-gluon or qq) have comparable energies. This convention 
deviates slightly from that of our previous publication [18| for which the scaling variable 
was the recoil mass of the gg and qq systems opposite the hard photon (M recoi i, defined 

by M recoi i = ^/4i^ eam (l — i? 7 /i?beam))- Comparison with continuum data taken ~20 MeV 
below each of the T resonances mitigates the effect of the ~1 GeV continuum center-of- 
mass energy extrapolation between the T(1S) and below-4S data samples required in the 
previous analysis[9j], for which continuum data were only taken in the 10.55 GeV center-of- 
mass region. To compare ggg with qq hadronization, we simply bin by scaled momentum of 
the particle in question. 



A. Event Selection 

We impose event-selection requirements identical to those used in our previous study 



of inclusive direct photon production in T decays [19f] . Those cuts are designed primarily 
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to suppress backgrounds such as two-photon collisions, QED events (including tau pair 
production), and beam-gas and beam-wall collisions. Luminosity, event count, and photon 
yields (z 7 >0.5) are given in Table fl] 

TABLE I: Summary of data and JETSET Monte Carlo used in analysis. For each data set, we 
track the number of photons per unit luminosity, as well as the total number of observed hadronic 
events per unit luminosity C HadEvts denotes the total number of events in each sample identified 
as hadronic by our event selection requirements. The number of photons having scaled momentum 
z 7 greater than 0.5 is presented in the last column. For BB Monte Carlo simulations, the small 
number of observed high-energy photons is a result of detector resolution and mis-reconstruction. 



Data Type 


Type 


Resonance 


E cm (GeV) 


C (pb- 1 ) 


HadEvts (xlO 3 ) 


N^(z > 0.5) (xlO 2 ) 


IS 


Data 


T(IS) 


9.455-9.465 


1220 


22780 


2190 


2S 


Data 


T(2S) 


10.018-10.028 


1070 


9450 


888 


3S 


Data 


T(3S) 


10.350-10.360 


1420 


8890 


795 


4S 


Data 


T(4S) 


10.575-10.585 


5520 


18970 


1650 


1S-CO 


Data 


< T(1S) 


9.400-9.454 


144 


515 


57 


2S-CO 


Data 


< T(2S) 


9.523-10.017 


312 


932 


103 


3S-CO 


Data 


< T(3S) 


10.083-10.349 


185 


532 


59 


4S-CO 


Data 


< T(4S) 


10.410-10.574 


2100 


5680 


647 


IS 


JETSET MC 


T(1S) 


9.455-9.465 




1160 


99 


2S 


JETSET MC 


T(2S) 


10.018-10.028 




9190 


700 


3S 


JETSET MC 


T(3S) 


10.350-10.360 




3890 


270 


4S 


BB MC 


T(4S) 


10.575-10.585 




8350 


3 


1S-CO 


JETSET MC 


< T(1S) 


9.400-9.454 




8170 


681 


2S-CO 


JETSET MC 


< T(2S) 


9.523-10.017 




7610 


666 


3S-CO 


JETSET MC 


< T(3S) 


10.083-10.349 




12850 


1150 


4S-CO 


JETSET MC 


< T(4S) 


10.410-10.574 




63630 


5680 



B. Background Suppression 

To determine the characteristics of resonant T — > ggj events, we must subtract the 
background arising from non-resonant qq'y and e + e~ — > tt^ events produced in continuum 
e + e~ annihilations at = Mr( n s)> with n=l, 2, or 3. This is done by direct scaling of the 
event samples collected off-resonance on the nearby continuum. 

In order to isolate continuum qtpf events, tt7 contamination must be explicitly sub- 
tracted, using a Monte Carlo simulation of tau pair events. We find that TT7 events com- 
prise about 5% of the qq^ data sample passing the event selection cuts[19|. Beam-gas and 
two-photon backgrounds were investigated and found to be negligibly small. The photon- 
tagged sample can also be contaminated by cases where the high-energy photon candidate 
is not produced directly, but is actually either a secondary daughter (mostly from 7r° decay) 
or a mis-identified hadronic shower. Figure [T] illustrates the fraction of photons in Monte 
Carlo simulations of on-lS resonance and below-4S continuum, respectively, that are not 
produced in a direct decay. 

Integrated over all tag photon momenta considered in this analysis, tt° contamination 
comprises a ~15% background to the direct-photon sample. Monte Carlo simulations also 
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FIG. 1: JETSET Monte Carlo prediction for fraction of photons not produced directly, but through 
the decay of neutral particles (such as ir°, rj, r/', and w) on the IS resonance (left) and on the 
continuum below the 4S resonance (right). 



indicate that the 7r° contamination tends to cancel when we take ratios of resonant photon 
production to continuum photon production. 



1. Particle Identification 

Our photon and particle identification procedures are identical to those developed in [l9| . 
Photon candidates are selected from showers with widths and patterns of energy deposition 
consistent with those of a photon, as opposed to a neutral hadron (e.g., tt° with overlapping 
photon showers, K^, neutrons, etc.). To ensure that the events are well-contained within 
the CLEO detector, we require |cos# 7 | < 0.707 (6> 7 defined as before as the polar angle 
between the beam axis and the direct photon). For p (and p), we require that charged 
tracks have specific ionizatation (dE/dx) and also RICH information consistent with those 
expected for protons. For momenta less than 1 GeV/c, we also require that the associated 
charged track dE/dx information be inconsistent (at the two standard- deviation level, with 
a the momentum-dependent specific ionization resolution) with that expected for true pions. 
Although this results in a discontinuity in particle identification efficiency at 1 GeV/c, this 
requirement is necessary to ensure a high-purity sample. For all p and p candidates, we 
require that p (p) momenta exceed 400 MeV/c to suppress beam-wall and fake backgrounds 
(i.e. K + and tt + that pass p identification cuts) and also to eliminate concerns regarding 
protons ranging out in the beampipe. For reconstruction of 4> (^2(1270)) from kaons (pions) 
we require that pairs of opposite charged tracks with momenta greater than 200 MeV/c (500 
MeV/c) have particle identification information consistent with their assumed identities. A's 
are identified using the standard CLEO algorithms for reconstruction of detached vertices. 
Tables HT1 and II III summarize the raw, observed particle yields for our measurements, for data 
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and Monte Carlo simulations, respectively. 



Particle Type 


(ggg) /(qq) [Data] 


(mi)/ (mi) [Data] 


A 

P 
P 

4> 

/ 2 (1270) 


(873600 ± 1400)/(107300 ± 600) 
(1399800 ± 1200)/(295900 ± 500) 
(1359500 ± 1200)/(285400 ± 500) 
(227900 ± 1600)/(48300 ± 800) 
(193000 ± 4000)/(66500 ± 1800) 


(3480±90)/(570 ±60) 
(7970±90)/(2190 ±50) 
(7830±90)/(2090 ±50) 
(1950± 150)/(380 ±70) 
(1600 ±400)/(400 ± 200) 



TABLE II: Data particle yields for the on-lS resonance compared to continuum events. First 
column is particle type. Second and third colunmns show particle counts for the data with in the 
format of (resonance yield) / (continuum yield) for the three gluon (2nd column) and two gluon 
one photon (3rd column) analyses. 



Particle Type 


(ggg)/(<E) [JETSET MC] 


(mi)/(<i<il) [JETSET MC] 


A 

P 
P 

<t> 


(136700 ± 500)/(1333200 ± 2000) 
(266600 ± 500)/(3334200 ± 1800) 
(257300 ± 500)/(3198300 ± 1800) 
(48100 ± 900)/(837000 ± 4000) 


(690±30)/(6410 ±150) 
(1650 ±40)/(20660 ± 140) 
(1590 ±40)/(19880 ± 140) 

(380±80)/(6000 ±800) 



TABLE III: Monte Carlo particle yields for the on-lS resonance compared to continuum events. 
First column is particle type. Second and third colunmns show particle counts for the data with 
in the format of (resonance yield) / (continuum yield) for the three gluon (2nd column) and two 
gluon one photon (3rd column) analyses. 



2. Backgrounds to the Proton Sample 

We use Monte Carlo simulations to assess fake proton backgrounds. Figure [2] illustrates 
proton fakes for a sample of below-lS Monte Carlo continuum simulations. The solid black 
curve shows the number of all particles identified as protons that were also tagged as true 
protons. The red dashed (blue dotted, magenta dash-dot) curve corresponds to those parti- 
cles that were identified as protons, but that were generated as true kaons (pions, positrons) 
in the Monte Carlo simulated event sample. Proton backgrounds are observed to be present 
at the ~10% level and are expected to largely cancel in the enhancement ratio. 

C. Signal Definition 

In this analysis we measure particle enhancements in both the ggg and gg^ decays of 
the Upsilon system, relative to (7(7(7) production on the underlying continuum. Our defi- 
nition of enhancement is given quantitatively as the continuum-subtracted resonance yield 
relative to the continuum yield. Thus defined, an enhancement of 1 indicates that a given 
particle is produced as often (per event) on the continuum as on the resonance. Note that 
our definition of 'continuum' here includes both continuum below the resonance peak, as 
well as resonance— >■ qq through vacuum-polarization; i.e., all e + e~ —>■ gg-like processes which 
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FIG. 2: (Left) Proton fakes for a sample of below-lS Monte Carlo simulations. The solid black 
curve shows the number of all particles identified as protons that were also tagged as true protons. 
The red dashed (blue dotted, magenta dash-dot) curve corresponds to those particles that were 
identified as protons, but that were actually kaons (pions, positrons). (Right) Same for on-lS event 
simulations. Note the discontinuity at 1 GeV/c, resulting from our momentum-dependent particle 
identification requirements below and above that momentum (see text). 

must be explicitly subtracted in determining the characteristics of 3-gluon resonant decays. 2 
Furthermore, note that for the T(2S) and T(3S) data, there is no subtraction of cascades 
to lower T states or Xb decays. In what follows, "T(2S)" denotes a sum over T(2S) direct, 
T(2S) — > T(1S)+X and Y(2S)^ jXb- Assuming the direct decays of the T resonances are 
identical, an T(2S) enhancement smaller than that of the T(1S) implies that the enhance- 
ments from the first and third processes enumerated above are therefore smaller than for 
the T(1S). 

In general we have two continuum-subtraction options: we may determine enhancements 
for all resonances relative to the below-4S continuum (for which the statistics are largest, 
but the difference in e + e~ collision energies is also largest) or we may find enhancements 
relative to their individual below-resonant continuua. For mass-fitted particles we normalize 
exclusively to the below-4S contiuum, as the individual continuua (below- IS, -2S, and -3S) 
have insufficient statistics to yield well-fitted mass peaks. For particle counts determined 
by the momentum spectra (protons and antiprotons), we normalize to both the below-4S 
continuum as well as the resonance-specific continuua and incorporate the differences in the 
enhancements calculated in the two cases into the overall systematic error. 



2 Vacuum polarization processes are subtracted by direct scaling of the continuum using the T — ► 7* 
values tabulated previously[19|. 
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D. Particle production in three-gluon vs. qq events 

The previous CLE0-I[9| analysis already observed significant enhancements of p and A 
produced in 3-gluon decays of the IS relative to the below-4S continuum. We repeat that 
analysis with our larger, current data set, as detailed below. Errors on particle yields are 
obtained from the error returned from the fit if the particle count is obtained by fitting a mass 
peak (A, 0, / 2 ), or by the square root of the total count if the particle count is obtained from 
a simple integration over the momentum spectrum (p, p). For the ggg analysis described 
below, we determine enhancements as a function of scaled momentum and also calculate 
momentum-integrated enhancements for each particle, to allow comparison with previous 
results. 



E. ggj Analysis 

For the ggj analysis we normalize the total particle yield to the photon count in a 
given photon momentum bin. For each bin, we then find the fractional contamination T 
of resonance photons "i?" due to the underlying continuum "C" (Eqn. [2]) in terms of the 
visible cross-section a for high-energy photons and the known beam energies U E" . 

Once T is known, the resonance yield can be extracted by straight-forward algebra. 



III. RESULTS FROM UPSILON DECAYS 
A. ggg Enhancements with respect to qq 

1. Baryon Enhancements 

Figure [3] presents our A enhancements binned according to scaled momentum, defined 
as before as the momentum of the particle divided by the beam energy. In the figure, 
blue square (gold triangle, green diamond) symbols correspond to enhancements on the IS 
(2S, 3S) resonance. Closed symbols are data and open symbols are JETSET 7.3j20[ event 
generator simulations followed by the full CLEOIII GEANT-based[2l|] Monte Carlo detector 
simulation. From the figure we see that the T's show qualitatively the same behavior for 
all resonances (IS, 2S, 3S) in both data and Monte Carlo, namely a smooth decrease in 
enhancement with increasing scaled momentum. We note that the enhancements decrease 
steadily as one goes from T(1S) to T(2S) to T(3S) and that the data, at all scaled momenta, 
show significantly greater enhancements than do the Monte Carlo simulations. 

Figure H] shows the p and p enhancements. With the exception of the very lowest mo- 
mentum bin, which is most subject to range-out effects, the consistency between the two 
indicates that beam-wall and beam-gas backgrounds (which produce an excess of p in the 
beam) are not substantial. As compared to A enhancements, p and p enhancements are 
lower and the differences between IS, 2S, and 3S enhancements (as well as the differences 
between data and Monte Carlo) are smaller. 
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FIG. 3: Raw (i.e., observed, and with no relative efficiency corrections applied) enhancements for 
ggg — > k + X binned according to scaled momentum (p\/ '-Ebeam)- Blue square (gold triangle, green 
diamond) symbols correspond to enhancements on the IS (2S, 3S) resonance. Closed symbols are 
data, open symbols are derived from JETSET 7.3 Monte Carlo simulations. No relative efficiency 
corrections have been applied to these 'raw' data. 
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FIG. 4: (Left) Raw enhancements for ggg — > p + X binned according to scaled momentum. Blue 
square (gold triangle, green diamond) symbols correspond to enhancements on the IS (2S, 3S) 
resonance. Closed symbols are data, open symbols are JETSET Monte Carlo. (Right) Same for p. 
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(gold triangle, green diamond) symbols correspond to enhancements on the IS (2S, 3S) resonance. 
Closed symbols are data, open symbols are JETSET Monte Carlo. (Right) Enhancements for 
/ 2 (1270). 



2. (j) and j \(1210) Enhancements 

Figure [5] shows <fi enhancement results binned according to scaled momentum. Symbols 
are as above with blue square (gold triangle, green diamond) corresponding to enhancements 
on the IS (2S, 3S) resonance. Closed symbols are data and open symbols are JETSET Monte 
Carlo. Here, we have normalized the resonant production at 9.46 GeV to the continuum 
production at 10.55 GeV. For production, the lowest momentum bins for the resonance 
are particularly sensitive to low-momentum kaon acceptance. Figure [5] also shows the f 2 
enhancement results binned according to scaled momentum. The f 2 peak is not well-defined 
at low momentum (lowest two bins). No Monte Carlo comparison is presented since our 
current Monte Carlo event generator, by default, will not generate / 2 tensor particles. 



3. Particle Momentum- Integrated Enhancements 

Figure [6] shows the particle enhancements integrated over all momenta for each particle, 
summarized numerically in Table W\ . We note that the baryons (A, p, p) have enhancements 
greater than 1, the meson enhancement is closer to unity, and the production of the tensor 
/ 2 is less than unity over our kinematic acceptance region. Our results are, in general, 
numerically consistent with the prior CLEO-I analysis, albeit with considerably improved 
statistical precision. 
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FIG. 6: Compilation of momentum-integrated enhancements for ggg events. Blue square (gold tri- 
angle, green diamond) symbols correspond to enhancements on the IS (2S, 3S) resonance. Closed 
symbols are data, open symbols are JETSET Monte Carlo. Systematic errors and relative effi- 
ciencies have now been included for this compilation. The CLE084 study did not measure an 
enhancement for /2(1270) and also only presented a single enhancement for the sum of p and p. 



B. gg7 Enhancements with respect to qq'j 

There are sufficient CLEO III statistics to present enhancements binned according to 
photon momentum, but integrated over particle momenta for A, p and p. For all particles, 
we also present momentum-integrated enhancements. 



1. Baryon Enhancements 

Figure [7] shows A results binned according to scaled photon momentum. For A's, as 
compared to the momentum- integrated ggg/qq enhancements, we observe a lower overall 
enhancement, on the order of 2 as opposed to 2.5-3 for the T's (Figure E]). Figure [S] shows 
p and p enhancement results binned according to scaled photon momentum. We note that 
ggi/QQl P an d P exhibit behavior similar to that of A's. 
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ments on the IS (2S, 3S) resonance. Closed symbols are data, open symbols are JETSET Monte 
Carlo. 
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FIG. 8: (Left) Raw enhancements for gg^ — > p + X binned according to scaled photon energy, 
integrated over all p momenta. Blue square (gold triangle, green diamond) symbols correspond to 
enhancements on the IS (2S, 3S) resonance. Closed symbols are data, open symbols are JETSET 
Monte Carlo. (Right) Same for p. 
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FIG. 9: Compilation of photon momentum-integrated enhancements for ggj events. Blue square 
(gold triangle, green diamond) symbols correspond to enhancements on the IS (2S, 3S) resonance. 
Closed symbols are data, open symbols are JETSET Monte Carlo. Systematic errors and relative 
efficiencies have now been included for this compilation. 

2. Photon Momentum- Integrated Enhancements 

Figure M shows the photon momentum-integrated enhancements for each particle, sum- 
marized numerically in Table |V] We note that all baryons show enhancements lower than 
in the 3-gluon case (Figure [6]). 



IV. INCLUSIVE PROTON PRODUCTION IN X bJ DECAYS 

Photon transitions of the T(2S) and T(3S) to the Xb( r ) states allow us to measure the 
baryon yields in Xfe(') decay, in association with a radiative transition photon 'tag'. Typ- 
ical photon tag energies in this case are of order 80-160 MeV. Due to the large 7r° — > 77 
backgrounds to such transition photons at these relatively low photon energies, which com- 
promises the statistical power of such tags, the data permit only an extraction of the proton 
and antiproton enhancements. Of particular interest is the proton yield in Xb2 vs. Xbi decays; 
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the former is expected to be dominated by decays via two gluons, the latter is expected to be 
dominated by decays to qq(g), with the gluon expected to carry away very little momentum. 

To ensure that photon-finding systematics largely cancel in the ratio, and to also exclude 
possible contributions from initial state radiation, we compare particle yields within the 
Xbj system directly rather than normalizing, e.g. relative to the underlying e + e~ — > 99(7) 
continuum. We first conduct a Monte Carlo study to determine the relative efficiency of 
reconstructing a J=2 transition photon relative to J=l event, and also the efficiency when 
we require that a proton be found in addition to the transition photon. We compile statistics 
on the XbO ~ > p(p) + analyses, separately for J=0/J=1 and for J=2/J=l. For the latter, 
the overlap of the two observed photon signals results in a highly correlated event yield 
for the two transitions. We correspondingly extract this ratio from a signal fit to a double 
Gaussian plus a smooth background. For the former, we simply fit two separate signal 
Gaussians directly. We find that the efficiency for reconstructing photon-proton correlations 
in Xb2 —> 99 decays is approximately 95% that for photon-proton correlations in Xbi — > <7<z(g) 
events. 

To check the sensitivity to our particle identification criteria, we have compared results 
using very tight proton identification requirements (with a reduction in efficiency by more 
than 50%) vs. the 'standard' loose proton identification criteria used above. We obtain a 
comparable correction factor for the J=2/J=l event yields using more restrictive particle 
identification criteria. 

Results are presented in Table IIVI We note that the observed enhancements are, again, 
smaller than those observed in comparing three-gluon fragmentation from the T resonance 
with qq fragmentation. 

V. CROSS-CHECKS AND SYSTEMATICS 

In order to verify our procedures and probe possible systematic uncertainties, two pri- 
mary cross-checks were employed. We first compare the Monte Carlo enhancements at the 
event generator-level with those determined after the generated events are processed through 
the full CLEO-III detector simulation ("detector-level"), as a function of momentum. In 
general, these enhancements will differ for several reasons, including differences in: a) the ef- 
ficiencies for finding recoil particles in qlpi vs. ggj events resulting from angular distribution, 
event multiplicity, and particle momentum differences, b) event selection efficiencies, c) tt° 
contamination levels, and d) recoil center-of-mass discrepancies between the continuum data 
under the T(1S) resonance vs. the below-T(4S) continuum. In cases where the generator- 
level and detector-level enhancements are statistically inconsistent with each other at the 
2<T level, we use the ratio between the generator-level and detector-level enhancements as a 
correction factor and take half of the amount by which this correction deviates from unity as 
an estimated systematic error. (Note that these corrections have already been incorporated 
into the results presented in Figures and [6] and [9} . Figures [10] and [11] shows the comparison 
of p enhancements determined at the event-generator vs. post-detector-simulation levels of 
Monte Carlo simulation. Integrated over momentum, typical corrections are typically of 
order 10%. 

In addition to the comparison of generator vs. detector-level enhancements, we have 
made an additional (largely redundant) check of possible biases due to non-direct photons 
resulting from, e.g., 7r° — > 77, 77 — > 77, etc. Using Monte Carlo simulations, we compare 
the enhancements obtained using direct-photons only, compared with the enhancements 
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FIG. 10: (Left) Scaled momentum binned enhancements for ggg — > p + X at generator level and 
after detector simulation. Blue square (gold triangle, green inverted triangle) symbols correspond 
to enhancements on the IS (2S, 3S) resonance. Closed (open) symbols are generator (detector) 
level Monte Carlo enhancements. (Right) Same for p. 



i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 

■ 1S Generator 

A 2S Generator 
|^ T 3S Generator 

□ 1 S Detector 

A 2S Detector 
V 3S Detector 



J I I I I I I I I I I I I I I I I I I I I I I L 



3031206-037 



~ i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



▼ 3S Generator 

□ 1 S Detector 

A 2S Detector 

V 3S Detector 



A 



_l I I I I I I I I I I I I I I I I I I I I I I L 



0.50 



0.60 



0.70 0.80 

X v 



0.90 



1.00 0.50 



0.60 



0.70 0.80 
X v 



0.90 



1.00 



FIG. 11: (Left) Scaled momentum binned enhancements for gg*y decays to p at generator level and 
after detector simulation. Blue square (gold triangle, green inverted triangle) symbols correspond 
to enhancements on the IS (2S, 3S) resonance. Closed (open) symbols are generator (detector) 
level Monte Carlo enhancements. (Right) Same for p. 
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TABLE IV: Summary of inclusive proton (and antiproton) results for Xbj decays. For checks 
of internal consistency, data have been separated into sub-samples, labeled with capital Roman 
letters. For J=2 relative to J=l, e.g., the scale of systematic uncertainties is set by the constancy 
of the value across sub-samples collected in different running periods (r.m.s.~0.03), the magnitude 
of relative efficiency corrections (~0.05) and the consistency of results obtained using different 
particle identification criteria. For summed results (labeled "all"), the second error shown is the 
systematic error. 
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tight 
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1.17±0.22 


(2S D) 


tight 


1.024+0.035 


0.88±0.20 


2S, all 




1.068 ±0.010 ±0.040 


1.11 ±0.15 ±0.20 


Monte Carlo (3S A) 


loose 


1.057 ±0.016 


1.030 ±0.072 


Monte Carlo (3S A) 


tight 


1.034 ±0.015 


1.042 ± 0.066 


Monte Carlo (3S B) 


tight 


1.041 ±0.013 


1.051 ±0.049 


MC, 3S all sets 




1.043 ± 0.008 


1.043 ± 0.036 


Monte Carlo (2S A) 


tight 


1.052±0.014 


1.121±0.058 


Monte Carlo (2S A) 


loose 


1.043±0.015 


1.076±0.061 


MC, 2S all sets 




1.046 ± 0.010 


1.061 ±0.025 



obtained when we include all Monte Carlo photons which pass our photon selection, inde- 
pendent of parentage. Integrated over momentum, this again constitutes a ~5% effect, and 
is conservatively included as an additional (in quadrature) systematic error. 

To test the sensitivity of our analysis procedures across different running periods, 
we have calculated the enhancements for photon-tagged T(4S) on-resonance events vs. 
photon-tagged below-T(4S) continuum events, spanning the full CLEO-III data set. Since 
Y(4S)^BB ~ 100%, we expect that any event having a photon with z 7 > 0.5 is a continuum 
event. Hence, the calculated enhancement should be zero. In all cases, save for p, we find 
good agreeement between the below-4S continuum particle yields per photon tag, and the 
on-4S particle yields per photon tag. For p, we find deviations from the null expectation 
at the level of ~5-7%, and incorporate these deviations (bin-by-bin in momentum) into our 
total systematic error for that particular case. For the case of the broad / 2 resonance, sen- 
sitivity to our parameterization of the smooth background also contributes a non-negligible 
systematic uncertainty. 

We note that most systematic errors cancel in our ratios. To summarize the sources of 
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systematic uncertainties, the largest components are the efficiency bias (~10%), the non- 
direct photon background (~5%), and the run- dependence of our result, as measured by 
our 'null' check (~6%), all added in quadrature. Our results with statistical and systematic 
errors are listed in Table EJ The statistical uncertainties in the data are typically of order 
10%, with the exception of ggj/qq^ for and f 2 , which are of poorer statistical quality. 

VI. DISCUSSION AND SUMMARY 

We have, for the first time, measured the momentum-dependent ratio of baryon and 
meson production in gluon vs. quark fragmentation at -^s ~10 GeV. After reproducing 
the previously measured per-event baryon production rates in three-gluon decays of the 
T(1S) resonance relative to the underlying continuum, we have extended that study to 
include the other narrow T resonances and, with higher statistics, now explicitly examine the 
momentum dependence of the enhancements for all these states. Integrated over momentum, 
we observe approximately 5% (10%) lower baryon production per-event for T(2S) (T(3S)) 
decays compared to the vector ground state. Nevertheless, the per-event production of A's 
for each of the narrow T resonances is observed to be greater than twice that of continuum 
fragmentation at the same center-of-mass energy. 

We additionally compare, for the first time, particle production in two-gluon vs. quark- 
antiquark fragmentation. We find, in particular, that baryon production (per event) in two- 
gluon decays is somewhat smaller (~20% for baryons) than that observed in three-gluon 
decays. For A production, we still observe a significant (~ x2) enhancement in two-gluon 
fragmentation relative to quark-antiquark fragmentation, although the excess enhancement 
for p is <10%. For p, which represent our highest-statistics sample, our results are incon- 
sistent with a model where baryon production in gluon fragmentation is only a function 
of the available center-of-mass energy; clearly, the number of fragmenting partons is also 
important, although our measured enhancements fall short of the expectations from a naive 
independent fragmentation model. Our results, for all measured integrated enhancements 
are presented in Table IVl 

Although event generators such as JETSET have had tremendous success in describing 
the gross details of particle production in e + e~ collisions, our study indicates that there may 
still be considerable tuning needed at the single-particle yield level. 
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Particle 


991 1 mi data 


991 1 mi MC 


999 /0Q data 


999 /qq MC 


A (IS) 
A (2S) 
A (3S) 


1.86 ±0.25 ±0.03 
1.98 ±0.27 ±0.08 
2.18 ±0.36 ±0.02 


1.38 ± 0.039 
1.38 ±0.018 
1.49 ± 0.023 


2.668 ±0.027 ±0.051 
2.333 ±0.019 ±0.021 
2.128 ±0.021 ±0.010 


1.440 ±0.003 
1.428 ±0.002 
1.450 ±0.002 


p(lS) 
P (2S) 
P (3S) 


1.21 ±0.11 ±0.03 
1.26 ±0.11 ±0.06 
1.51 ±0.17 ±0.06 


1.582 ±0.034 
1.495 ±0.018 
1.53 ±0.021 


1.623 ±0.014 ±0.116 
1.469 ±0.011 ±0.103 
1.348 ±0.013 ±0.116 


1.331 ±0.005 
1.177 ±0.003 
1.214 ±0.003 


p(lS) 
p(2S) 
P(3S) 


1.45 ±0.14 ±0.26 

1.46 ±0.12 ±0.17 
1.39 ±0.17 ±0.27 


1.589 ±0.034 
1.513 ±0.018 
1.51 ±0.020 


1.634 ±0.014 ±0.111 
1.500 ±0.011 ±0.102 
1.323 ±0.013 ±0.115 


1.333 ±0.005 
1.175 ±0.003 
1.210 ±0.003 


4> (is) 

0(2S) 
0(3S) 


1.78 ±0.49 ±0.08 
1.73 ±0.52 ±0.06 
1.87 ±0.81 ±0.06 


0.673 ±0.013 
0.658 ±0.012 
0.662 ±0.015 


1.423 ±0.051 ±0.065 
1.308 ±0.041 ±0.041 
1.355 ±0.054 ±0.047 


0.836 ± 0.003 
0.805 ± 0.001 
0.808 ± 0.002 


/ 2 (1270) (IS) 
/ 2 (1270) (2S) 
/ 2 (1270) (3S) 


1.34 ±0.84 ±0.15 (< 2.74) 
2.22 ±1.53 ±0.20 (< 4.68) 
1.41 ± 1.48 ±0.10 (< 3.87) 




0.658 ±0.058 ±0.175 
0.621 ±0.094 ±0.171 
0.702 ±0.104 ±0.175 





TABLE V: Numerical summary of momentum-integrated enhancement results. Second and third 
columns show results from Figure El the photon momentum integrated ggi/qqi study. Fourth and 
fifth columns show results from Figure [6j the particle momentum integrated 999 /qq study. MC 
refers to JETSET Monte Carlo. Data errors are statistical and systematic; MC errors are purely 
statistical. For the / 2 , we present 90% C.L. upper limits, given the poor statistical significance of 
the ggi/qqj enhancements. 
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